Strong field effects on pulsar arrival times: circular orbits and equatorial beams 
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ABSTRACT 



If a pulsar orbits a supermassive black hole, the timing of pulses that pass close to the hole 
will show a variety of strong field effects. To compute the intensity and timing of pulses that 
have passed close to a nonrotating black hole we introduce here a simple formalism based on two 
"universal functions," one for the bending of photon trajectories and the other for the photon 
travel time on these trajectories. We apply this simple formalism to the case of a pulsar in circular 
orbit that beams its pulses into the orbital plane. In addition to the "primary" pulses that reach 
the receiver by a more-or-less direct path, we find that there are secondary and higher order 
pulses. These are usually much dimmer than the primary pulses, but they can be of comparable 
or even greater intensity if they are emitted when pulsar is on the side of the hole furthest from 
the receiver. We show that there is a phase relationship of the primary and secondary pulses that 
is a probe of the strongly curved spacetime geometry. Analogs of these phenomena are expected 
in more general configurations, in which a pulsar in orbit around a hole emits pulses that are not 
confined to the orbital plane. 



Introduction 



Pulsars are rotating neutron stars that emit beams of radiation that can be detected on Earth. As 
the star rotates, this beam sweeps past the Earth, producing regular pulses. The timing of these pulses 
is tied to the large inertial moment of a compact body, making it an extreme ly stable clock: puls ars have 
been found whose pulse arrival times fluctuate by less than 200 nano-seconds (jVerbiest et al.l 120081 ). Pulsar 
timing is thus an excellent probe of delicate phenomena, including gravity, in the vicinity of the pulsar. It 
has been used to detect the tug of planets orbiting a pulsar, to observe the gradual loss of orbital energy to 
gravitational waves in binary sys tems containing a pulsar, and has been proposed as a meth od of detecting 
cosmological gravitational waves ( Estabrook fc Wahlquist 19751 : Sazhin 19781 Detweiler 1979h . In this paper 
we explore how pulsar timing might probe the gravitational environment of pulsars orbiting supermassive 
black holes. 

Supermassive black holes are now thought to be ubiquitous in the Universe, residing in the cores of 
most large galaxies, including our own. They have masses in the ra nge of millions to billions of Solar masses : 
our own Galaxy' s black hole has an estimated mass of 4 x 10 6 M Q ( Ghez et al. 2005 : Eisenhauer et al. 20051 
Ghez et al.l 120081 ) . The best estimates of its mass come from observing the orbits of O-type stars in the 
Galactic nucleus. These observations have also revealed that the Galactic nucleus is home to a significant 



-2 - 



population of young massive stars, contrary to earlier expectations (|Lu et al.ll2006T ). Perhaps the best model 



for this population has it forming in situ from the dense molecular hydrog en disk surrounding the black hole 



with an initial mass funct ion that is strongly tilted towards massive stars (jNavakshin fc Sunvaevll2005t iLevi 



20071 : iManess et al.ll2007l) . The Galactic nucleus is thus a likely environment for neutron stars to form, some 
of which may be pulsars orbiting deep within the potential well of the supermassive black hole. Although 
no such pulsars have been detected, they may be discovered by future searches for periodic signals with 
corrections for high accelerations (due to gravity) and/or dispersion (due to plasma in the Galactic core). 

The systems of interest here will consist of neutron stars at distances down to within several million 
kilometers (a few Schwarzschild radii) of the central black hole. Even at these close separations, it would take 
many years for the neutron star orbit to decay due to gravitational radiation emission, making it possible to 
conduct lengthy timing observations of a pulsar deep within the strong field of a black hole. (By contrast, a 
pulsar orbiting within a thousand Schwarzschild radii of a IOMq black hole would decay within a year.) 

In this paper we consider a pulsar to be in orbit around a supermassive black hole, and we ask what 
effects the strong field of the hole would have on pulsar observations. Such effects would depend on the 
details of the hole/pulsar system, and there are many details: the pulsar orbital elements, the alignment 
of the pulsar spin axis and the orbital plane, the pulsar spin rate, the angle between the pulsar spin axis 
and the pulse emission direction, and the black hole spin. (The black hole mass can be treated as a scaling 
parameter.) 

An important step in understanding the details of strong field effects on pulsar obsrvations is to un- 
derstand them in simple cases. This is what we do here, in two stages. First, we consider pulses emitted 
from a pulsar in the neighborhood of a Schwarzschild (nonrotating) black hole, with no restrictions on the 
pulsar orbit, spin, or emission direction. We point out that in the case of a spherically symmetric hole it 
is not necessary to do extensive computing of null geodesies. Rather, it is necessary only to compute two 
functions of the emission direction, one representing the bending of the light path, and the second the time 
delay along that path. These two "universal functions" of emission direction are parameterized only by the 
distance of the pulsar from the hole at the emission event. 

The second stage in our analysis is to apply these universal functions to a particularly simple astro- 
physical scenario, that of a pulsar emitting in the orbital plane as it travels in a circular orbit around a 
Schwarzschild hole. It turns out that even in this simplest possible case, there is a rich set of interesting, and 
potentially observable phenomena. Analogs of these phenomena could occur also in binary pulsar systems 
and in binary neutron star/black hole systems of comparable mass. 

The paper is organized as follows. In Sec. [2] we define the two universal functions, and we present 
quantitative results for them. Next, in Sec. [3K we show how we can use these functions to find any timing 
effect in the Schwarzschild geometry. We then, in Sec. [3Jd, apply this method explicitly to pulses beamed 
in or near the orbital plane. These results are then applied, in Sec. [4J to the simple case of a circular orbit 
and equatorial beaming. In particular, we show effects that would appear to a distant receiver of the pulse 
trains. A summary and a consideration of future applications of this method arc given in Sec. [5] To allow the 
main ideas of the paper to be as clear as possible, several sets of details have been relegated to appendices. 
Appendix A gives the details of finding the direction of emission as a function of time for a pulsar in a 
relativistic circular orbit. Appendix B gives the details of the computation of the universal f unctions. The 



relativistic calculations presented use the notational conventions of the text by Misner et al. (jMisner et al 



19731 ): in particular we use c = G = 1. 
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2. The universal functions for deflection and timing 




Fig. 1. — Photon path from radial coordinate ry to R. The path on the right indicates mild bending of the 
path due to spacetime curvature. On the left is shown the highly bent photon path which has an impact 
parameter near the critical value of photon capture, and a 4> ou t, at large R, of 2tt. 



In a Schwarzschild hole with spacetime metric 

ds 2 = - (1 - 2M/r) dt 2 + (1 - 2M/r)~ 1 dr 2 + r 2 ( d9 2 + sin 2 6 d(/) 2 ) , (1) 

we consider a point at r = To, (f> = 0, in the 9 — ir/2 equatorial plane. At this point, a photon (i.e. , a 
null geodesic) with 4-momcntum p' 1 is emitted in the equatorial plane. From the fact that p$ and pa are 
constants of motion for the equatorial photon, it follows that the photon orbit obeys 

1 / dr\ 2 1-2M/V 1 



,,4 



where b is the photon impact parameter. 

We define 0; n to be tan _1 (rp^/p r ), the angle with respect to the outgoing direction, at which the photon 
is emitted. (Note that this differs from the angle tan -1 [rp^ / ' yj 'g rr p r ) that would be measured in the local 
frame of a coordinate stationary observer.) The value 0i n is a constant characterizing the orbit, and is related 
to the impact parameter b by Eq. @ with r set equal to tq: 



tan 2 0; n r b 2 



A critical value of the impact parameter is & cl -it = 3v3 M; for b < fe C rit, inward going photons will be 
captured by the hole. For any r , capture will occur if |0j n | > cr it where C rit is the root, between ir/2 and 

7T, of 

( r 2 2MY 1/2 
Wc rit = -(^-l + — J • (4) 

For \(f>i n \ less than cr it, our universal functions relate the description of the photon at very large distance 
R to its emission conditions. The first of these functions measures the bending of a photon trajectory, relative 
to the radially outgoing direction from its emission point. To define this function we start by defining 
J-(4>i n ;R) to be the <f> location, at large radius i?, of a photon emitted at radius ro, at angle 0i n to the 
outgoing direction. We then define 

F(0 in ) = lim T(4> in ;R) . (5) 

R — >oo 

The function F is an odd function of 4>- m ; it is parametrized by r /M, but has no other dependencies. The 
second universal function T is related to the coordinate time T(0; n ; R) required for the the photon to reach 
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asymptotically large R. To have the result be a finite value, we subtract the coordinate time to reach R with 
0in = 0, so that 

T(0 in ) = lim [T(0 in ; i?) - T(0; i?)] , (6) 

i? — >oo 

and is an even function of </>j n ; like F, it is parametrized by ro/M, but has no other dependencies. 



250 py 




Fig. 2. — Universal functions F and T. Curves are labeled with the values of ro/M: 5, 10, and 100. The 
values of </> crit , from Eq. ((J) are respectively 2.0895, 2.6109, and 3.0896. 

Plots of F and T are given in Fig. [2] for three values of ro/M. Approximate fits to these curves, good 
both for small <f>- m and for </>i n near ^> C rit, are given by 

F(M = -signOfe) log (l - \M) - + <j>- m (7) 

V 0crit / <Pcrit 

T(4> in ) = r„ (1 - cos(0 in )) - 3V31og (l - ^) - ^-\<f> in \ ■ (8) 

V 0crit / 9crit 

The details of the computation of these curves, and more accurate analytic approxmations to the curves, are 
given in Appendix B. 



3. Pulse arrival times from universal functions 
3.1. General case 

We consider an emission event, at Schwarzschild radial coordinate ro, and a photon emitted in the 
direction n. We let a be the angle, as shown in Fig. [3l between the orbital plane and the plane containing 
the photon trajectory, i.e. , the plane defined by n and the radially outward direction. (If these directions 
coincide, we take a to be zero.) 

We introduce a spatial triad of orthonormal vectors e x , e y ,e z , with e^ej, in the pulsar orbital plane, 
and we define e x > , e y > , e z < , to be the orthonormal spatial triad that result from a rotation by a around the 
outgoing direction e x , so that e X ',e y i lie in the plane of the photon trajectory. 

To specify the direction of n we use the spherical polar angles <f>o, 0o with respect to the e x ,e y ,e z 
directions. (Note that for convenience we are using the latitude j3, rather than the more typical colatitude 
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: 7r/2 — /3.) From the rotational transformation between the unprimed and primed triads we get 

tan 0o 



tan /% 

tana = - 

sin 4>q 

cos 4>- m i = cos /?o cos 4>o 



(9) 
(10) 
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Fig. 3. — Triads and angles used to define photon directions. 

Here 0j n < is the angle, in the plane of the photon trajectory, between the photon direction and the 
radially outgoing direction. We can, therefore, find t/w , the asymptotic (f> direction, from 

0oo' = F(&n0 ■ (11) 

spherical polar angles, in the e x ,e y ,e z frame, for the final photon direction, are 

cos — sin a sin c/)oa> tan^oo = cosatan^oo' . (12) 

We summarize here the set of equations that give the final photon direction /Jooj^ooj and relative time of 
arrival t 00l in terms of the original photon direction /3o, </>o : 

(13) 
(14) 

(15) 
(16) 



COS 4>in' — COS /?o COS 



sin Poo = sin 



>(*»o) , a sin/3 ; a 

V sin po + cos 2 po sin ^ 

tan^oo = tan (F((/>in0) i cos /3 sm ^ 

V sin /3 + cos 2 /3 sin </) 

ioo = *e + T(0 in /) . 



3.2. Special case: beaming in orbital plane 

If the pulsar beam is emitted into the orbital plane, then po = and the relations in Eqs. (fT3|) - (fT6|) 
reduce to <fri n / — <po, and (f>oo = F(4>i n ). The photon remains in the orbital plane, so Poo = Pq = 0. 

For the simplest example of strong field effects we restrict ourselves to the case in which the receiving 
antenna is precisely in the orbital plane of the pulsar. This means that photons reaching the antenna must 
have exceedingly small values of Pq (of order of the antenna size divided by tens of kiloparsecs). In this case 
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the initial-to-final transformations in Eqs. (|13p -(|16 p can be replaced by approximations valid to first order 
in P : 

sm 0o 

0oo = Ffa). (18) 

In these equations it is assumed that /3q -C 4>q, a condition that will be violated only by photons directed very 
nearly radially outward, and hence experiencing no strong field effects. The sine function in the numerator 
of Eq. (jTTJ) implies /3oo = if F(0o) = i". This corresponds to the case that all small-/?o photons with the 
same value of tfio are focused into the orbital plane, creating an intensity amplification that would be infinite 
aside from diffraction limitation. 

Equations (fTT)) . (fTS)) give us the factors by which strong field effects lead to convergence or divergence 
of the photon beam both in the pulsar orbital plane / p i an c, and perpendicular to it / pcrp : 

, _dF_ _ sin TO )) ( . 

/plane — , , Jpcrp — . , ■ 

a<po sm tpo 

The total strong-field effect on the cross sectional area of the beam reaching the receiver is / p i an e/perp- The 
effect on the photon flux at the receiver will therefore be 1/ (/piane/pcrp)- The total radio intensity reaching 
the receiver will also depend on the gravitational redshift of the photons and the Doppler shift due to the 
motion of the emitting pulsar. Both effects are of order M/tq. 



4. Appearance and timing of pulses 

We make another simplifying assumption: that the pulsar that emits its beam in the equatorial plane, 
is moving in a circular orbit at radius ro. For definiteness we take the angular location of the pulsar to be 
given by O rb = $lt, where "i" here and below, indicates coordinate time. At a particular emission time t e , 
it is shown in Appendix [A] that the angle at which the photon is emitted is given by 

r O (1 - 2M/r )- 1/2 + cos (7^ - cj]t + S) 

tan (</> in (t c ) = , — — — ■ — . 20 

v ' 7 1 sm(7 L [il — uj\t c + 0) 



Here ro is the radial coordinate of the circular orbit; = y/M/r^ is the pulsar orbital angular velocity (per 
unit coordinate time); u> is the pulsar spin rate as measured by a comoving observer; the Lorentz factor 7 is 
- 3M/r ; and S is a phase constant specifying the direction of the beam at t c = 0. 

If a distant radio receiver is at re c, then to find the emission time of photons that are destined to be 
received, we must solve 

nt e + F((P in (t e )) =0 rec . (21) 

Due to the nature of F near the critical angles (f> ct it, there are, in principle an infinite number of solutions 
corresponding to photon orbits that circle the gravitating center zero times, once, twice, etc. It will be useful 
to refer to these received pulses as primary, secondary (once around the gravitating center), tertiary, and so 
forth. We shall see, however, that these distinctions can become ambiguous. 

The observationally important questions are the timing and appearance of the pulses that show strong 
field effects. We start by arguing that the effect on the pulse shape will be negligible if the rotation period 
(seconds or less) of the pulsar is much smaller than the pulsar orbital period (1000s of seconds to years, 
depending on r /M and on the mass of the central supermassive black hole). 
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Fig. 4. — The pulsar beam sweeping through the narrow range of directions that connect to the receiver. 



Figure [4] shows the geometry of beam emission, with Ar/>b cam indicating the inherent angular width (in 
the orbital plane) of the pulsar beam, and with A</()p r i mary and A0 sccon d ary indicating (but exaggerating) the 
receiver acceptance directions, i.e., the range of photon directions that connect to the distant radio receiver. 
The angular size of A(/> pr j mary is approximately the ratio of the receiver diameter to the receiver distance. 
The much smaller angular size A0 secon d ar y is further divided by dF/d(f>i n ^ a large number. For tertiary and 
subsequent beams, the same description applies except that the value of dF/d4>in is even larger. 

If the pulsar were at a fixed coordinate position, and rotating at w, then the shape of the received pulse 
can be viewed as the result of the narrow cones of receiver acceptance sweeping through the A^bcam beam 
profile. This viewpoint makes it clear that the time profile of the primary, secondary, etc. beams would 
have the same shape. The pulsar, of course, is not coordinate stationary, but is orbiting with orbital speed 
Q. This means that there will be a small change in the nature of dF/d<j>i n during the passage of the beam 
width through the receiver acceptance cones, since the value of <t>- m for received photons will change slightly 
during beam reception. To estimate this effect we can consider the change AF during the change in pulsar 
orbital location A0 ol b: 

dF . dF n Fo SI 

— Ar> in « — A0 orb - . Pa e F 2tt - . (22) 

dtpin d<p- m IjJ LO 

(The second approximations assumes that 0; n « Gr j t and that F is large.) For a secondary beam, F must 
be of order 27r; for a tertiary beam, F must be order 47r etc. Thus, in the case of high order beams or for 
exceptionally large values of Q/ui, there could be some distortion of received pulse shapes. This possibility 
will not be considered further in the current paper. 

Figure [5] shows a schematic of pulses emitted from a pulsar orbiting at a distance 30 x GM/c 2 from 
a supermassive black hole of mass M. The pulsar rotation rate has been chosen to be very low, only 20.2 
times the orbital frequency, in order to show clearly some of the phenomenology of the pulse arrival times. 
For a 4 x 10 6 M Q black hole, this means an orbital period of 2tt x 30 3//2 x GM/c 3 ps 20,000 seconds and 
a pulsar rotation period around 1000 seconds; an actual observed pulsar would likely be rotating hundreds 
or thousands of times faster. The horizontal axis represents the pulse arrival time (relative to an arbitrary 
start time), while the vertical axis gives the attenuation factor 1 //piano due to horizontal spreading of the 
beam. Using this factor as a tag on the pulses is useful in the discussion of the pulse sequences. The factor 
gives a rough indication of relative pulse strengths, though the full amplification or attenuation of the beam 
requires the complete spreading factor l// p i an o/pcrp, and will be discussed below. 

Initially, the taller set of pulses represents pulses that arrive at the detector from the pulsar along a 
more-or-less direct path. The initially weaker but strengthening set represents pulses that arrive from a path 
that is bent around the black hole in a sense prograde to the orbit. Figure [S] illustrates how this secondary 
path "unwinds" and becomes more direct as the pulsar moves along its orbit. The effective photon path 
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Fig. 5. — The factor l/f p \ ane for photon capture events from a pulsar with orbital radius 30 x GM/c 2 , and 
pulsar rotation rate 20.2 times the orbital rotation rate. Events for several orbits are shown on the left. 
Details are shown on the right for the orbital epoch during which the pulsar is on the far side of the hole. 



length £ also shortens with time, giving rise to a shorter (blueshifted) pulse period P = P (l + £/c). Around 
time 600 x GM/c 3 , the pulsar passes behind the black hole (the middle panel in Fig. [5]). After this, the 
prograde path becomes the more direct path, while what was formerly the more direct path now gets wound 
around the black hole in the retrograde sense, causing its pulses to weaken and the period of the pulses to 
redshift. The pulse timing near the time of this orbital phase is shown in greater detail in the second panel 
of Figure where the prograde- wound pulses are shown in gray to distinguish them from the other pulses. 
The process repeats itself one orbit later. The ragged pulses at the bottom of the two graphs are pulses from 
paths with higher-order windings around the black hole; some of these will eventually unwind to become, 
for a time, the most direct path. 

Figure [7] shows the same system, but with the vertical axis now showing the full pulse amplification 
(or attenuation) factor 1 / /plane /perp- Around the same time as the "primary" and "secondary" pulse trains 
swap roles, both sets of pulses go through a spike of amplification due to strong lensing around the black 
hole. 

It should be understood that the "intensity" in Figure [7] and below refers to the photons received per 
unit time. The true energy intensity would include the effect of Doppler shifts due to the orbital motion. 
These effects, of order r$Q,/c, are significant - around 20% for an orbit with tq = 30 x GM/c 2 - but are 
omitted to emphasize the photon path effects. 

The details of the pulse timing and amplification phenomenology depend on the orbital radius. To 
illustrate this, Fig. [5] shows the same information as Fig. [5] but for a more relativistic system, where the 
pulsar orbital radius is only 10 x GM/c 2 . Again, the pulsar rotation rate is taken to be artificially slow 
(only 47 times the orbital frequency) in order for the figure to show distinct pulses. In this more relativistic 
system, the pulse period shift, beaming asymmetry, and "tertiary" (multiply- wound) pulse paths are more 
apparent. 

Figure [5] shows the pulse trains that a radio astronomer might actually see from the tq = 30 x GM/c 2 



to observer 



Fig. 6. — Prograde and retrograde photon trajectories. The pulsar is shown orbiting in the counterclockwise 
direction around a much more massive black hole. The leftmost cartoon shows that the observing radio tele- 
scope receives a primary (P) direct pulse and a highly bent secondary (S) pulse. The subsequent panels show 
trajectories as the pulsar continues its orbtial motion. The primary trajectory has increased bending and 
the secondary less, until the pulsar is directly opposite the receiver and the two trajectories are symmetric. 
The following panels shows how the prograde trajectory then becomes the primary (less bent) trajectory. 




UN 

) 900 



Time of Arrival (GM/c 2 ) Time of Arrival (GM/c 2 ) 



Fig. 7. — The relative intensity of pulses for the same system shown in Fig. [5] 

system of Figs. [5] and [71 but assuming a more realistic pulsar rotation rate that is thousands of times faster 
than the orbital frequency. For such a sufficiently fast rotation rate, the pulsar orbital position changes 
negligibly during the emission of a set of a few pulses. The pulse intensity, determined by F(^>i n ), and 
dF/dipin, and the red/blue shift in the period, determined by £, therefore are fixed, once the orbital phase 
is fixed. This means that we can show pulse intensity and primary vs. secondary phase shifts without 
specifying any one pulsar rotation rate. For this reason, Fig. [5] does not specify units of the horizontal axis; 
it is understood that the spacing between direct pulses is the period of pulsar rotation. The time that is 
specified for each panel is the time that determines the orbital position when the set of pulses is emitted. 

The vertical axis represents the pulse power relative to an undeflected beam, and each pulse is given 
a Gaussian shape. Panel (a) is a segment of the pulse train when the pulsar is on the near side of the 
black hole, and only a single train of direct pulses is visible. Panel (b) is a segment from the orbital phase 
during which the pulsar is passing behind the black hole: a secondary train of prograde-bent pulses, with 
blueshifted period, starts to appear. In panel (c), the prograde path has become the more direct path, 
and the pulses from what was originally the direct path are weakening, and their period redshifting, as the 
path becomes a bent retrograde one around the black hole (though both sets of pulses are amplified due to 
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Fig. 8. — The same plot as in Fig. [5l but for a more relativistic system, with with orbital radius 10 x GM/c 2 , 
and with rotation rate 47 times the orbital frequency. 

lensing). In panel (d), the original pulse train has disappeared almost entirely, leaving only the pulses from 
the newly-unwound path. 

We note that when the most-direct path shifts from passing retrograde around the hole to prograde 
around the hole (see Fig. [6]), there is a sudden jump in the frequency of the strongest pulse train, as the 
derivative of the effective path length t switches from lengthening to shortening. (Equivalently, in this 
geometry, a train of high-frequency "secondary" pulses rises up and replaces the low-frequency "primary" 
pulses.) This is the strong-field analog of the well-known cusp in the Shapiro time delay curve, where the 
time delay switches from increasing to decreasing; since the observed pulse period is multiplied by 1 plus 
the time derivative of the time delay P h s = -Po(l + dTshnpim/dt), this corresponds to a discontinuous jump 
in period. Unlike the case of the "standard" Shapiro delay, the discontinuity is finite even for a perfect 
alignment, since the deflecting mass is a black hole of finite size. There is also a slight asymmetry due to 
the special relativistic beaming ("headlight effect") of the pulsar in its orbit. 

5. Summary and Conclusions 

We have introduced the problem of black hole effects on pulses from an orbiting pulsar. We have shown 
that in the case of a spherically symmetric (i.e., nonrotating) hole, the use of two "universal functions" 
removes the need for extensive computation of null geodesies. The universal function approach to compu- 
tation has been used for a first exploration of the phenomena that might be encountered with pulses from 
a pulsar-hole system. This first exploration investigates pulse emission in the orbital plane, for a pulsar in 
circular orbit around a nonrotating hole. 

Even for this highly simplified configuration we have found a number of interesting phenomena: (i) "Pri- 
mary" pulses (pulses that travel from the pulsar to the receiver with relatively little gravitational bending) 
are continually accompanied by higher order (secondary, tertiary,...) pulses emitted during earlier pulsar 
orbits, (ii) As would be expected, the primary pulses do not have a constant period, but rather have a 
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Time of Arrival (arbitrary units) near t r = 100 GM/c 3 
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(b) Time of Arrival (arbitrary units) near t r = 540 GM/c 3, 
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Time of Arrival (arbitrary units) near t r = 600 GM/c 3 
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Time of Arrival (arbitrary units) near t r = 1000 GM/c 3 



Fig. 9. — The appearance of received radio pulses for several orbital epochs of a pulsar orbiting with radius 
30 x GM/c 2 . The pulsar rotation rate is assumed to be much greater than the orbital frequency. 



period that is modulated by orbital motion, (iii) The period of primary and of higher order pulses is not 
precisely the same, thus the higher order pulses arrive with a phase shift, relative to the primary pulses, that 
varies from one pulse to the next, (iv) Strong field effects dominate the pulsar observations when the pulsar 
is on the far side of the hole, the side opposite that of the receiver. In this case the emitted pulse can reach 
the receiver only by a highly bent path, (v) During the epoch of emission from the far side, the sequence of 
primary pulses and the sequence of secondary pulses, exchange roles. One consequence is that there is no 
clear distinction of primary and secondary for pulses emitted when the pulsar is close to the middle of its 
passage through the far side of the hole, (vi) For emission during most of the pulsar orbit, the secondary 
pulses have much lower intensity than the primary pulses. As the pulsar moves toward the dark side, how- 
ever, the secondary pulses increase in intensity and the primary pulses decrease, (vii) Among the pulses 
emitted from the dark side are pulses that are amplified by strong field effects analogous to gravitational 
lensing. 

The simple configuration studied would be expected to exaggerate strong field effects when compared to 
a more realistic configuration in which pulse emission is directed well out of the orbital plane. For example, 
if the pulsar spin axis is perpendicular to the orbital plane (as in our simple configuration), and emission 
is not close to the orbital plane (contrary to our simple configuration), the pulse trajectory will never cross 
the orbital plane, and hence never come sufficiently near the hole for strong field effects to be significant. 
More interesting is the case in which the spin axis and beaming details are such that the beam does cross 
(or pass close to) the orbital plane. Particularly noteworthy would be a configuration in which the receiver 
receives no primary pulses, and receives pulses only with the aid of strong field bending. It is also of note 
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that some of the phenomena we have described, especially the existence of higher order pulses, can also occur 
in principle in binary pulsar systems, and pulsar-hole binaries of comparable mass. 

Work is underway on investigating configuratins with out of orbit beaming. An exploration of the large 
parameter space will be feasible with the efficiency provided by the universal function approach, so we are, 
at least at first, restriciting attention to nonrotating supermassive holes. 
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A. Emission direction as a function of time 




Fig. 10. — The figure illustrates the orbital motion of the pulsar at frequency Cl, the parallel transported 
spatial direction , and the pulsar beam direction rotating at frequency lo relative to V , as observed 
in a frame comoving with the pulsar. 



In the orbital plane we let the angular position of the pulsar be <fi — Qt, with t the Scharzschild coordinate 
time and the pulsar's 4-velocity components are 

U° = 7 = 7« (Al) 

where 

7 - 1 = 1 . (A2) 

y/1 - 2M/r - r 2 W y/l - 3M/r 

We let V be a 4- vector that is spatial (that is, orthogonal to U^), that has no component out of the 
equatorial orbital plane, and that is parallel transported around with the pulsar. (Since the pulsar world 
line is a geodesic, this is the same as Fermi- Walker transporting V\) It is straightforward to show that the 
components of such an "inertial direction" are 

V r =j- 1 r Ksm(- f - 1 nt + S) V* = K cos {j^flt + S) V° = ° K cos fr^Qt + 5) . (A3) 

where K is a scaling constant and <5 is a phase factor determining the direction in which is pointing at 
t = 0. We let be a spatial vector, in the equatorial plane, that points in the direction of the pulsar 
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beam. We let ip be the angle, measured in the positive sense (the positive sense for fi) from V* to . In 
terms of the proper time r measured by the pulsar, we define the locally observed pulsar spin rate lu by 

ip = ut = uit/'y . (A4) 

From W^U^ = and cos?/> = WV^/VW we get 

A A 

W r = — sin (7- 1 [n - w]t + 5) = — cos (7- 1 [VL - u]t + 8) 

7 r 

W ° = l-2M/r C ° S [Q " ]t + S) ■ (A5) 

In the comoving frame of the pulsar the spatial direction of the photon beam is W*, thus the photon 
4-momentum must have the form p = kU + W. The value of k follows from p^p^ — and we find 

p= — ; A U + W, (A6) 

7v /l - 2M/r 

from which we get the components 

p° = - = + - r °"f ; cos - u]t + 5) (A7) 

y/1 - 2M/r l-2M/r K ' 

1 An a , 

y/l - 2M/r r Q 

p r = j^Aami-f-^n-^t + S). (A9) 
This tells us that the photon starts out its journey with 

p r dr roj sin(7 [O — u>]t + 6) 



with 



and with 



d ^ +cos r 1 -i[Q_ LJ ]t + 8) ' 

y/l-2M/r W 1 



p r _ dr _ 7- 1 sin(7" 1 [fi - u]t + 5) 

~f~~di~ , 1 + . :°" cos(7- 1 [»-^]t + (5) 

^/l-2M/r l-2M/r J ' 

pV d(f> ^/l-2M/r r « W L 1 ' 



(A10) 



(All) 



Equation (|20|) follows from Eq. (|A10[) and from the definition of tan0j n as the initial value of rdcp/dr. 



(A12) 



B. Computation of the universal functions 
B.l. Integral for tp^ 

For a photon beamed in the equatorial plane, we assume that we know its initial radial location r$ and 
its initial direction, as specified by the angle measured with respect to the outgoing radial direction 

(f) in = tan" 1 (r d<p/dr) . (Bl) 
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The equation for the photon orbit is 



1 fdr\ 2 l-2M/r 1 

r 4 I dd> r 2 b 2 ' K ' 



From the known initial values of r and of dr/dip we solve Eq. (|B2|) for the impact parameter b, a constant 
of motion for the orbit. We proceed with the following steps to find </>oo. In the following we describe 
as "delicate" any integral that has a divergent integrand, and for which special techniques, desribed in 
Appendix IB. 31 must be used. 

1. Suppose cos (0in) > 0, so that the photon is emitted going generally outward. In this case the photon 
will be moving only to larger r, but it may be moving to larger or smaller 0, and we must separate 
this case into two subcases depending on the initial value of dr/d<f>: 

(a) If initially dr/d<f) > 0, then we use 

f°° dr 

(j) 00 (t)=m+ ; (B3) 

y ' Jro ^r i /b{t) 2 -r 2 + 2Mr V ' 

(b) If initially dr/d<j> < we use 

f°° dr 

</•«,(*) = Ot- / ; (B4) 

V ; Jr ^Jr i /b{t) 2 -r 2 + 2Mr V ' 

In both these cases, the integrands do not diverge, and no special techniques are needed to carry out 
the integration. 

2. If cos0j n — it means that the photon is emitted on a trajectory tangent to the pulsar orbit. In this 
case we must check whether the the photon is going in the direction of increasing or decreasing <f>: 



(a) If sin</>i n > 0, then we use Eq. (|B3[) . 

(b) If sin^in < 0, then we use Eq. ([Elf . 

In the cos<?!>i n = case, either of the integrals is "delicate" since the integrand diverges at r = ro- 

3. If cos(0; n ) < and if b 2 > 27M 2 , then there will be a value of r, less than ro, at which dr/dcf) = 0, 
i.e., at which the denominator of the integrand in Eq. (|B3j) vanishes. Call that value r m i n . Then we 
need to consider the usual subcases: 

(a) If sin ((/>j n ) > 0, then 

f ro dr r°° dr 

(j )oo {t)=ttt + 2 + / - (B5) 
Jrmin ^Jr A /b{t) 2 -r 2 + 2Mr J r „ ^Jr±/b{t) 2 - r 2 + 2Mr 

(b) If sin(0 in ) < 0, then 

C r ° dr f°° dr 

4 >00 (t)=ttt-2 - / - (B6) 
7r mta y/ryb(t) 2 -r 2 +2Mr J ra ^Jr±/b{t) 2 - r 2 + 2Mr 

4. If cos (</>;„) < and if b 2 < 27M, then the photon will be captured by the black hole and there is no 
meaning to 0oo . 
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B.2. Integral for time to infinity 

We need to find the coordinate time it takes for a photon to reach "infinity." The calculation must be 
divided into two subcases depending on whether the photon starts going outward (cos (</>in) > 0) or inward 
(cos (<pin) < 0). In both cases the calculation is based on the equation for dr/dt for the photon motion (see 
MTW, Eqs. 25.64 and 25.66): 




1. If the photon starts its trajectory at time t e , with cos (4>i n ) > 0, then the photon will always be traveling 
outward and 

roc 

too = U + / J(r)dr . (B8) 

J r 

In this case the integral is not delicate. 

2. If the photon starts its trajectory with cos (fan) = 0, then Eq. (|B8j) again applies, but now the integral 
is delicate, since /(ro) is infinite. 

3. If the photon starts its trajectory with cos (<fi- m ) < 0, then let r min have the same meaning as in Eq. (|B5|) 
(|B6[) . The time of arrival of the pulse is now given by 

pro poo 
+ 2 / J(r)dr + / J(r)dr . (B9) 

« ^min J ro 

The first integral is delicate since J(r m i n ) is infinite. 



B.3. Computational details for integrals 



The integrals for </>oo and too cannot in general be expressed as elementary functions, and must be 
handled by numerical methods. If there is no singularity in the integrand, such as in the case cos(0i n ) > 0, 
we use adaptive Simpson quadrature with an absolute error limit set at 10~ 6 . 

In the case that the integrand does diverge, the divergences occur at the a root of the polynomial 
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2Mb 2 . For b 2 < 27 M 2 this polynomial has a negative root and two postive roots. We denote the 



largest root as r m i n since it is the minimum radius the photon trajectory will reach. We then write 

r 4 /6 2 - r 2 + 2Mr = (r - r min )/(r) (B10) 

where / = r 2 (r + r m ; n )/6 2 — 2Mr /r m - m does not vanish in the intervals of integration that are considered. 
The use of this is illustrated for one of the integrals occurring in Eqs. (|B5[) and (|B6|) : 



dr 



lin v / r 4/ 6 2_ r 2 + 2Mr 
dr 



dr 



(Bll) 



1 V/( r min) 



dr 
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dr 



dr 



ma \/ r ~ r min\/ /(''mill) 

The numerator in the second integrand is 
f(r 

so that the evaluation of I is reduced to 



/(r mi „) - f(r) 



) - f(r) = - p(r - r min )(r 2 + 2r mln r + 2r min - 26 2 M/r min ) 



7 = 2 



1 



V r ~ r min 

y/ f(r min) b 2 ^f f{r min ) 



Vr - r min (r 2 + 2r min r + 2r 2 nin - 2b 2 A//r mill )rfr 

(/7(^) + v / 7M)v / 7W 



(B12) 



(B13) 



The remaining integral is nonsingular and can be evaluated e.g., with an adaptive Simpson's routine. 

An additional complication arises if b 2 is very close to 27M 2 . In this case the two positive roots of /(r) 
approach each other, and the integral becomes very large. To deal with this case we introduce the notation 
r 2 for the smaller positive root and r 3 for the negative root, and we write 

dr 




tin "\Z^*min (^*n 

dr 



r3)(r- r min )(r - r 2 ) 
dr 



yjr(r - r 3 )(r - r min )(r - r 2 ) \/r m in{r 



dr 



- r 3 ){r - r min )(r - r 2 ) 
V r = r min(r + r min - r 3 ) dr 



)(r r 2 ) 7r min ^y/ r min (r m in - r 3 ) + y/r(r - r 3 ) >/r(r - r 3 )(r - r 2 ) 



\J ^min( r min ~~ 



2 log 



'2 



r 3 ) dr 



' r »«" (\/r m i„(r m i n - r 3 ) + ^r(r - rjjn ^Jr(r - r 3 )(r - r 2 ) 

(B14) 

The remaining integrand is well behaved at r — r m i n and is straightforward to evaluate numerically. 



To get asymptotic approximations for q 
3Af + and r 2 w 3Af - yT/3. If we set c, 

(HI) we find that A</> ex e and hence cx (r n 



near <^ cr ; t we take b 2 = 27M 2 + e, and we find that r min w 
= </>crit — A0, and assume positive ^ then from Eqs. j3]) and 
— T2) 2 . The logarithmic dependence in Eq. (|B14p becomes 

(B15) 



log(l/Vr m in - r 2 ) f» logA0 W -|log(l - 0in/0crit) ■ 

The prefactor of the logarithm in Eq. (IB14j) is evaluated to 2 when the approximations b = 3\/b, r m [ n = 
r 2 = 3M, r 3 = — 6M are used, and we get the asymptotic approximation — ^log(l — (f>i n / 4'ci-it) ■ The first 
term in Eq. {7|) follows from using this asymptotic approximation for Eq. (|B14j) in Eq. (|B5j) and modifying 
the result to make it an odd function of <p- m . The remaining terms are added to Eq. (Q to give F{<p- m ) ps <p- m 
for small </>;„. Equation (JT]), therefore, gives the right qualitative behavior of F(cj>i n ) both for </>;„ — » and 

0in — * </>crit- 

The approximation for T(0 m ) in Eq. j7]) follows from similar considerations with the addition of a 
Roemer time delay ro(l — cos </>i n ), the time it takes the pulsar signal (with c = 1) to cross the orbit on its 
way to the receiver. The appearance of this effect is clear in the t-q = 100A/ curve for T(<j>i n ) in Fig. [2] 



It has been useful to have more accurate approximations to the universal functions than the expressions 
in Eqs. ([7]) and ([8]). To that end we have fit the errors (differences from the numerical integrals) in Eqs. ([7|) 
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and (|8j with even-order Chebyshev polynomials. The following polynomials give better than 1% fits to the 
numerical calculations when used with the coefficients in Table [T] 



-F(0in) = -sign((/> in ) 



io g (i - Jpi) + f 2( f>l + Mi + Ml + Ml 



(B16) 



T(<kn) = r (l - cos0 in ) - 3V3 log(l - ^) - ^ |«fc n | + t 2 cj>l + U<f>i + t 6 0f n + t s $ n + t w ^ . (B17) 
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Table 1. Coefficients of residual of universal function F and T 

rp(M) </>crit h h fa h U *e *s tip 

5 2.0895 0.0634022 -0.0707841 0.0160038 -0.000411158 -0.86385 0.930642 -0.160082 0.00547166 0.0 

10 2.6109 0.0171645 0.00324181 -0.000180664 0.000315048 -0.263976 0.599057 -0.115977 0.00609665 0.0 

30 2.9677 -0.0352832 0.0569467 -0.0110258 0.000910931 -0.414603 0.786804 -0.212809 0.0236697 -0.00100843 

100 3.0896 -0.128706 0.134169 -0.0274787 0.00192051 -0.884133 0.929001 -0.23247 0.0249892 -0.00103073 



